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Abstract: Present data do not rule out a light neutral Higgs particle with mass 
\ below 40-50 GeV in the framework of 2HDM with tan/3 ~ 20-30. The promising 

possibility of searching for a light Higgs particle in such a scenario in photoproduc- 
tion at HERA collider is discussed. For the MSSM there is only a very small chance 
to observe the Higgs sector, and only for limited mass range ~ 45-50 GeV and with 
^ ! large tan (3. 

o\ : 
o . 

1 Introduction 

^ : 

The possibilities of Higgs searches at the HERA collider have been studied in the first HERA 
Workshop in 1987 [||]. It was found that it 'is (almost) impossible' to observe the Standard 
,-£! • Model (SM) Higgs. A similar conclusion for the SM scalar Higgs search at HERA, even with an 
upgraded luminosity and/or proton energy, can be found in the contribution to this Workshop 



X 



■ Non-minimal Higgs boson production at HERA has also been investigated during the HERA 



Workshop'87 as well as in other papers 113, 14, 15]. It was found that photon-gluon fusion into 



bbh(A) may be an important production mechanism of Higgs bosons in the two Higgs doublet 
extension of the SM at HERA. Also other subprocesses with Higgs boson bremsstrahlung, 
namely those with the resolved photon in the initial state, are important at HERA []TJJ]. Another 
production mechanism is gluon-gluon fusion via a quark loop, where the Higgs particle is 
produced in resonance. For large tan f3 and Higgs masses below 30 GeV this process dominates 



the production cross section over fusion ITS 



According to LEP I data, Higgs bosons in the MSSM have to be heavier than 45 GeV, 
therefore their production rate at the HERA collider is rather small. For the mass of h and 
A equal to 45 GeV both the jg and the gg cross sections for tan/3=30 are ~ 5 fb. When 
adding the similar contribution from WW fusion into h , one may expect 20-30 events to be 
produced at HERA with an integrated luminosity of 1 fb -1 . 

The situation is quite different in the non-supersymmetric version of the two Higgs doublet 
extension of the SM (the so called general two Higgs doublet model - 2HDM) since in this model 
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one light neutral Higgs boson with mass below 40-50 GeV may still exist, moreover with large 
coupling to r and b -quark. In this case there is a good chance to study the Higgs sector at 
HERA, with one thousand events expected for a Higgs mass of 5 GeV and tan/?=30, assuming 
£ ep =250 pb- 1 . 

Below we present the status of this model, i.e. , the 2HDM with a light Higgs boson, 
and in the next section we discuss the possibility to perform a Higgs boson search at HERA, 
focusing mainly on the gluon-gluon fusion production via a quark loop. Existing limits from 
LEP I for the coupling of a light neutral Higgs boson in the 2HDM are rather weak. Also 
present data on (g — 2) for the muon improve only slightly the limits on tan/3 for a Higgs 
mass below 2 GeV. Therefore it is extremely important to check if more stringent limits can 
be obtained from HERA measurements. The combined exclusion plot showing the potential of 
HERA measurement is presented in Sec. 4. Sec. 5 contains our conclusion. 

2 Status of the 2HDM with a Light Neutral Higgs 

The mechanism of spontaneous symmetry breaking proposed as the source of mass for the 
gauge and fermion fields in the SM leads to a neutral scalar particle, the minimal Higgs boson. 
According to the LEP I data, based on the Bjorken process e + e~ — > HZ*, it should be heavier 
than 66 GeV @. A possible extension of the SM is to include a second Higgs doublet to the 
symmetry breaking mechanism. In the two Higgs doublet models the observed Higgs sector is 
enlarged to five scalars: two neutral Higgs scalars (with masses Mh and Mh for the heavier and 
lighter particle, respectively), one neutral pseudoscalar (Ma), and a pair of charged Higgses 
(Mh±). The neutral Higgs scalar couplings to quarks, charged leptons and gauge bosons are 
modified with respect to analogous couplings in SM by factors that depend on additional 
parameters: tan/?, which is the ratio of the vacuum expectation values of the Higgs doublets 
V2/V1, and the mixing angle in the neutral Higgs sector a. Also further new couplings appear, 
e.g. Zh(H)A and ZH + H~ . 

In the following we will focus on the appealing version of the models with two doublets 
("Model II") where one Higgs doublet with vacuum expectation value t> 2 couples only to the 
"up" components of fermion doublets while the other one couples only to the "down" compo- 
nents fTj . (In particular, fermions couple to the pseudoscalar A with a strength proportional to 
(tan f3) ±l whereas the coupling of the fermions to the scalar h goes as ±(sina/ cos/3) , where 
the sign ± corresponds to the isospin =Fl/2 components). In such a model FCNC processes are 
absent and the p parameter retains its SM value at the tree level. Note that in such a scenario 
the large ratio v%/vi ~ m top /mb ^> 1 is naturally expected. 

The well known supersymmetric model (MSSM) belongs to this class. In the MSSM there 
are additional relations among the parameters required by super symmetry, leaving only two 
free parameters (at the tree level) e.g. Ma and tan/3. In the general case, denoted 2HDM, 
masses and parameters a and (3 are not constrained. Therefore the same experimental data 
may lead to very distinct consequences depending on which version of the two Higgs doublet 
extension of the SM, supersymmetric or nonsupersymmetric, is considered. 
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2.1 Present constraints on the 2HDM from LEP I. 



Important constraints on the parameters of the two Higgs doublet extensions of the SM were 
obtained in the precision measurements at LEP I. The current mass limit on the charged 
Higgs boson M H ±=A4 GeV was obtained at LEP I from the process Z — > H + H~, which is 
independent of the parameters a and (3. (Note that in the MSSM version one expects M#± > 
My/)- For the neutral Higgs particles h and A there are two main and complementary sources of 
information at LEP I. One is the Bjorken process Z — > Z*h which constrains g\zz ~ sin 2 (a— 
for M h below 50-60 GeV. The second process is Z — > hA, constraining g 2 ZhA ~ cos 2 (a — 0) for 
Mh + Ma < Mz- Results on sin 2 (a — (3) and cos 2 (a — (3) can be translated into limits on the 
neutral Higgs bosons masses Mh and Ma- In the MSSM, due to th additional relations among 
the parameters, the above data allow to draw limits for the masses of individual particles: 
M h > 45 GeV for any tan/5 and M A > 45 GeV for tan/? >1 @, §. In the general 2HDM the 
implications are quite different, here only the large portion of the (M^M^) plane, where both 
masses are in the range between and ~50 GeV, is excluded M. 

The third basic process to search for a neutral Higgs particle at LEP I is the Yukawa 
process, i.e. the bremsstrahlung production of a neutral Higgs boson h(A) from a heavy fermion: 
e + e~ — > ffh(A), where / means here b quark or r lepton |22|, || A new analysis of the 
Yukawa process by the ALEPH collaboration Jl0| led to an exclusion plot on tan/? versus 
the pseudoscalar mass, Ma- (The analysis by the L3 collaboration is also in progress [PH.). 
However, the obtained limits are rather weak allowing for the existence of a light A with mass 
below 10 GeV with tan/3 = 20-30, for M A =40 GeV tan/5 till 100 is allowed! For the mass 
range above 10 GeV, similar exclusion limits should in principle hold also for a scalar h when 
replacing the coupling tan/3 — > sin a/ cos/3. However, one would expect larger differences in 
the lower mass region, where the production rate at the same value of coupling for the scalar is 
considerably larger than for the pseudoscalar. More stringent limits should be obtained there@. 

In the following we will study the 2HDM assuming that one light Higgs particle may exist. 
Moreover we will assume according to LEP I data the following mass relation between the 
lightest neutral Higgs particles: Mh + Ma ^ Mz- We specify the model further by choosing 
particular values for the parameters a and (3 within the present limits from LEP I. Since 
sin 2 (a — (3) was found ||, [| to be smaller than 0.1 for 10 < M h < 50 GeV, and even below 
0.01 for a lighter scalar, we simply take a = (3. This leads to equal strengths of the coupling 
of fermions to scalars and pseudoscalars. Note that then the EW gauge boson couplings to the 
Higgs scalar h disappear For the scenario with large tan/3 ~ 0{m t jm^) a large enhancement 
in the coupling of both h and A bosons to the down-type quarks and leptons is expected. 

Below we present how one obtains limits on the parameters of the 2HDM from current muon 
(g — 2) data P3| , also the potential of the forthcoming E821 experiment |23| is discussed. (See 
Ref.ra for details.) 



2.2 Constraints on the 2HDM from (g - 2), 



The present experimental limit on (g — 2) for the muon, averaged over the sign of the muon 
electric charge, is given by [EB| : 



a 1 — 
% - 2 



A does not couple to W and Z [fjj. 
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The quantity within parenthesis, a exp , refers to the uncertainty in the last digit. 



The theoretical prediction of the SM for this quantity consists of the QED, hadronic and 
EW contributions: 



,SM 



a QED + a had + a EW 



The recent SM calculations of a M are based on the QED results from [29(1, the hadronic con- 



tributions obtained in |33], [31], [J^, [36]] and J^J and the EW results from p7[ . The 
uncertainties in these contributions differ among themselves considerably (see below and in 
Ref. JZS, |28], [H], 2~0|). The main discrepancy is observed for the hadronic contribution, therefore 
we will here consider case A, based on Refs. [p9|, 3~0| , |3~4 , |33| , |36| , |28|| , with relatively small error 



in the hadronic part, 
case B (Refs. [0, |3T 



It corresponds to 



,SM. 



= 1 165 918.27 (0.76) -10" 9 . (The results for 



|) with the 2 times larger error in the hadronic part is discussed in 



The room for new physics, like the 2HDM with a light scalar or a light pseudoscalar, is basi- 
cally given by the difference between the experimental data and the theoretical SM prediction:^ 
a e * p — a^ M = 5a^. Below, the difference 5a M for the considered case, A, is presented together 
with the error cr, obtained by adding the experimental and theoretical errors in quadrature (in 
1(T 9 ): 

4.73(8.43) and lim±(95%) 



-13.46 < 8a„ < 19.94 



One can see that at the 1 a level the difference Sa^ can be positive or negative. For the beyond 
the SM scenario where the contribution of only one sign is physically accessible (i.e. positive 
or negative Sa^), the other sign being unphysical, the 95%C.L. limits should be calculated [|(| 
separately for the positive and for the negative contributions (lim±(95%) above ). 

The future accuracy of the (g — 2)^ experiments is expected to be a™™ ~ 0.4- 1CT 9 or better 



24] , One also expects an improvement in the calculation of the hadronic contribution 
such that the total uncertainty will be basically due to the experimental error. Below we will 
assume that the accessible range for the beyond SM contribution will be smaller by a factor of 
20 compared to the present (Zim±95%) bounds. So, we consider the following option for future 
measurements (in 10~ 9 ): 



lim ± neu, (95%) : -0.69 < 5a n ™ < 1.00. 



The difference Sa^f w we now ascribe to the 2HDM contribution, so we take 5a^ 
and 5a n a ew 
scenarios: 



a 



{2HDM) 



a^ HDM ^ for present and future (g 



2) M data, respectively. We will consider two 



a) pseudoscalar A is light 

b) scalar h is light. 



Here we calculate the 2HDM contribution assuming for case a) a^ HDM ^ (Ma 



'(Ma), 



whereas for h) : a^ HDM ^ (M^) = a^(M/j). This simple approach is based on the LEP I mass 
limits for charged and neutral Higgs particles and differs from the full 2HDM predictions, stud- 
ied in Ref. f20| , significantly for a Higgs mass above about 30 GeV. Note that the contribution 



for the scenario b) is positive, whereas for the scenario a) it is negative. 

3 However in the calculation of a^ w the (SM) Higgs scalar contribution is included (see discussion in [^0] 
4 An improvement in the ongoing experiments at low energy is expected as well. 



4 



The exclusion plots for tan/5 obtained from present (g — 2)^ data at 95%C.L. for a light h 
or A, go beyond those from LEP I for Higgs masses below 2 GeV. tan (3 ~ 15 is still allowed 
for the mass of the Higgs particle as low as 1 GeV, above 2 GeV tan (3 is limited to 20. These 
results together with others will be presented later in Sec. 4. 

It is worth pointing out the unique role of the forthcoming (g — 2)^ measurement in clar- 
ifying which scenario of the 2HDM is allowed: the model with a light scalar or with a light 
pseudoscalar. If the 5a^ ew is positive (negative) then the light pseudoscalar (scalar) is excluded. 
Further constraints on the coupling of the allowed light Higgs particle can be obtained from 
HERA, which is very well suited for this task. 



3 Search for a Higgs particle at HERA 

We now study the possibility of light neutral Higgs scalar and/or pseudoscalar production in a 
2HDM at HERA [14], [0| |16| . We limit ourselves to the mass range above 5 GeV, in order to 
apply the LO approach. Next order results are in preparation [[LB], but we expect that even a 
K-factor ~ 1.5-2 will not change the results drastically. The results obtained for 2HDM hold 
also for MSSM, provided the proper range of mass is considered, i.e. , above 45 GeV. The 
results relevant for SM can also be obtained from the 2HDM predictions for h production with 
tan (3=1. 

Photon-gluon and gluon-gluon fusions in photoproduction at HERA are expected to be 
basic sources of a light Higgs bosons in 2HDM |T5[]. Note that the ZZ and WW fusions are 
not relevant since the pseudoscalar, A, does not couple to EW gauge bosons, and the scalar 
couplings, which are proportional to sin(a — /3), are put to zero by the assumed by us condition: 
a = (3. (Even if the upper experimental limits for sin(a — f3) are used (Sec. 2.1), light Higgs 
boson production in ZZ /WW fusion is still effectively suppressed.) 

The total cross section for on-shell neutral Higgs boson production is calculated, along with 
the rates for the particular final states — > t + t~ or bb. These decay channels are the most 
important since in 2HDM with large tan/3 h and A decay mainly to the heaviest available 
fermionic "down- type" states. (Details can be found elsewhere, e.g. , in §, |T7|1). 



3.1 Bremsstrahlung of the Higgs boson: 
lated processes. 



7<7 — > bbh(A) and other re- 



At HERA, the photoproduction of neutral Higgs boson in a 2HDM via 

75 -> bbh(A) 



(1) 



may be substantial [ 13| , 15| . The total h cross section (integrated over the bb final state) is 
presented in Fig. la. Note that this process also includes 76 — > bh(A), as well the lowest order 
contributions due to the resolved photon, i.e., bb — > h(A), bg — > h(A)b, etc. These subprocesses 
were studied in Ref. [|14|, in Fig.2a we present obtained results. Each of these processes needs 



an independent analysis of the background [18]. As this work is not yet completed, we will not 
use these processes for the derivation of the exclusion plot in Sec. 4. 
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HERA ep (346 GeV):gg->h and gamma g->bbh 



HERA ep (314GeV):gg->h (A) 




M_h,GeV 



*s\ tan beta = 30 

\ h(A)->all 




M_h(A),GeV 



Figure 1: a) The total cross section for h production at HERA with the upgraded energy. The 
gg fusion (solid lines) and ^g fusion (dotted lines) are considered. The results for tan/5 =30 
(upper curves) and the tan/3=l (lower curves) are shown. GRV parametrizations were used 
for both the photon and the proton [18]. b) The cross sections for gg fusion for the nominal 
HERA energy (with tan P =30). Total h (solid upper line) and A (dashed upper line) production 
cross-sections are shown, along with the results for the t + t~ final state (solid and dashed lower 
lines for h and A, respectively). 



3.2 Gluon-gluon fusion via a quark loop 

A Higgs boson photoproduction in gluon-gluon fusion via a quark loop, 

99 - h(A), (2) 
can be even more significant ITl, |I5|. The results for HERA with upgraded energy are also 



presented in Fig. la. A comparison of (2) with 7$ fusion (1) shows that, for large tan/3 and 
for mass below ~ 30 GeV, the gg fusion clearly dominates the total cross section. Note that 
the total gg cross section for tan (3=30 is large: a ~ 10 5 fb for a Higgs mass of 5 GeV, falling 
to 5 fb at mass ~ 45 GeV (where the a gg meets the o~ ig ). 

A mass of 45 GeV corresponds to the lowest currently allowed mass for MSSM Higgs bosons. 
(Note that in MSSM, h and A tend to be degenerated in mass for large tan/?). Adding the 
contribution for tan/3=30 from processes (1) and (2) for both scalar and pseudoscalar as well 
as that due to WW fusion into h M which are of the same order, we estimate that HERA will 
produce 20-30 events of this sort with luminosity 1 fb -1 . 

In Fig. la the tan/3=l case corresponds to the prediction for SM Higgs production in 7$ 
and gg fusion. Applying the current limit for the SM Higgs mass, Mmggs > 66 GeV || , we see 
how small the corresponding 7^ and gg cross sections are, more than three orders of magnitude 
smaller than the rate found in WW fusion into h @. 

Comparing the scalar production in Fig. la with corresponding production in Fig. lb, in 
which the nominal HERA beam energy was used, we see an extremely weak dependence on 
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beam energy for scalar production. In Fig. lb we also consider pseudoscalar production via 
gluon-gluon fusion (2) for tan/3=30. The total cross sections (h(A) — > all) and the cross 
sections for the t + t~ final state are shown. It is interesting to notice the large difference, in the 
mass range below 10 GeV, between scalar and pseudoscalar production seen both in the total 
cross section a as well as in the (T-Br(— > t + t~). Note that the difference almost disappears for 
the r + r~ cross section above the bb threshold. 




Figure 2: In the ^yp CM system at y/s w =170 GeV a) the cross section for scalar production 
in various subprocesses (tan /3 =20, from Ref. j^j ) and b) the rapidity distribution for t + t~ 
pair are presented. The background (77 — ► t + t~ ) and the signal (the scalar Higgs boson 
contribution, integrated over AMh=l GeV) are shown (from Ref.^T3{J). 

For detection, it is useful to study the rapidity distribution da/dy of the Higgs bosons in the 
7P centre of mass system. Note that y = —\log = — \l°9^ri where x p {x 1 ) are the ratio 
of the energy of the gluon to the energy of the proton and photon, respectively. The (almost) 
symmetric shape of the rapidity distribution found for the signal is extremely useful to reduce 
the background and to separate the gg — > h(A) contribution, which we will discuss now for the 
t + t~ final state. 

The main background in the mass range between t + t~ and bb thresholds is due to 77 — > 
t + t~ . In the region of negative rapidity da/dy T + T - is very large, e.g. ,for jp energy equal 
to 170 GeV the cross section ~ 800 pb at the edge of phase space (y T + T - ~ —4), and it then 
falls rapidly when y T + T - approaches 0. At the same time, the signal reaches at most 10 pb (for 
Mh=5 GeV). The results for a scalar are shown in Fig.2b. The region of positive rapidity is 
not allowed kinematically for this process since here one photon interacts directly with x 7 = 1, 
and therefore y T + T - = —\log— < 0. A significantly different topology found for 77 — > t + t~ 
events than for the signal should allow a reduction of this background. 

The other sources of background are qq — > t + t~ processes (not shown here). These processes 
contribute to positive and negative rapidity y T + T -, with a flat and relatively low cross-section 
(below 0.5 pb) in the central region. 
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Note that Higgs decaying into fe-quarks has a much more severe background, and we will 
not discussed it here (see Ref.fLU p|). 

Assuming a luminosity £ ep =250 pb _1 /y we predict that gg fusion will produce around one 
thousand events per annum for = 5 GeV (and roughly 10 events for = 30 GeV). A 
clear signature for the tagged case with a t + t~ final state at positive centre-of-mass rapidities 
of the Higgs scalar should be seen, even for a Higgs mass above the bb threshold (more details 
can be found in Ref.[|l^]). For the pseudoscalar case even more events are expected in the mass 
region below 10 GeV. 

To show the potential of HERA, the exclusion plot based on gg fusion via a quark loop with 
the r + r~ final state can be drawn. In this case, as we mentioned above, it is easy to find the 
part of phase space where the background is hopefully negligible. To calculate the 95% C.L. 
for allowed value of tan (3 we take into account signal events corresponding only to the positive 
rapidity region (in the 7^ CM system). The results for the ep luminosity C ep =25 pb _1 and 
500 pb _1 are presented in Fig. 3 and will be discused in the following section. 



4 Exclusion plot for 2HDM 

In Fig. 3 the 95% C.L. exclusion curves for the tan/3 in the general 2HDM ("Model II") ob- 
tained by us for a light scalar (solid lines) and for a light pseudoscalar (dashed lines) are 
presented in mass range below 40 GeV. For comparison results from LEP I analysis presented 
recently by ALEPH collaboration for pseudoscalar |T(J is also shown (dotted line). The region 
of (tan/3, Mh(A)) above curves is excluded. 

Constraints on tan/5 were obtained from the existing {g — 2) M data, including LEP I mass 
limits (Sec. 2. 2). We see that already the present (g — 2) M data improve LEP I limits on tan/? 
for Ma ^ 2 GeV. A similar situation should hold for a 2HDM with a light scalar, although here 
the Yukawa process may be more restrictive for Mh < 10 GeV 0. 

The future improvement in the accuracy by a factor of 20 in the forthcoming (g — 2) M 
experiment may lead to more stringent limits than provided by LEP I up to a mass of h or A 
equal to 30 GeV, if the mass difference between scalar and pseudoscalar is ~ M z , or to even 



higher mass for a larger mass difference |20| . Note, however, that there is some arbitrariness in 
the deriving the expected bounds for the 5a™ ew . 

The search at HERA in the gluon-gluon fusion via a quark loop may lead to even more 
stringent limits (see Fig. 3) for the mass range 5-15 (5-25) GeV, provided the luminosity will 
reach 25 (500) pb _1 and the efficiency for the r + r~ final state will be high enough The other 
production mechanisms, such as ^g fusion and other subprocesses with the resolved photon, 
are expected to improve these limits further W\ IT5L |TS| . 



In the very low mass range, additional limits can be obtained from the low energy NL 



77 collider with v /s^=10 GeV. In Ref.fl7H we found that the exclusion based on 77 fusion 
into Higgs, decaying into n + fi~, may be very efficient. In Fig. 3 the results corresponding to a 
luminosity 10 fb _1 are presented. 



5 In this analysis the 100% efficiency has been assumed. If the efficiency will be 10 % the corresponding limits 
will be larger by factor 3.3 (a simple scaling) for fixed luminosity. 
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Exclusion plot (95 % C.L.) 
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Figure 3: The 95% C.L. exclusion plot for a light scalar(solid lines) or light pseudoscalar 
(dashed lines) in 2HDM. The limits derivable from present (g — 2) M measurements and from 
existing LEP I results (pseudoscalar production in the Yukawa process) (dotted line) are shown. 
The possible exclusions from HERA (the gluon-gluon fusion via a quark loop with the t + t~ 
final state) for luminosity 25 pb~ x and 500 pb' 1 as well from 77 — > at low energy NLC 

(10 fb -1 ) are also presented (from [17J). Possible limits from future data for (g — 2) M are also 
shown. The parameter space above the curves can be ruled out. 



5 Conclusion 



In the framework of 2HDM, a light neutral Higgs scalar or pseudoscalar in the mass range 
below 40-50 GeV is not ruled out by the present LEP I and (g — 2)^ data. The other low energy 
experiments cover only part of parameter space of 2HDM; some, such as the Wilczek process, 
have large theoretical uncertainties both due to the QCD and relativistic corrections ([12|, [7]]) (see 



also discussion in [1T5T |T6l 



The role of the forthcoming (g — 2) /t measurement seems to be crucial in clarifying which 
scenario of 2HDM is allowed: with light scalar or with a light pseudoscalar. If the ba n ^ w is 
positive(negative) then the light pseudoscalar(scalar) is excluded. Then, further constraints on 
the coupling of the allowed light Higgs particle can be obtained from HERA, which is very well 
suited for this. The simple estimation performed at luminosity 500 pb _1 for one particular 
production mechanism, namely gluon-gluon fusion, is already promising; adding more processes 
may further improve the situation significantly. The most important experimental handle is a 
good efficiency for the t + t~ channel. 

All this suggests the large discovery /exclusion potential of HERA for the mass range 5-20 



GeV ||15|| . It is unlikely that the LEP/LHC experiments will have a larger potential in such a 



mass region [[14] . 

The very low mass region may also be studied at low energy NLC machines. We found 
that the exclusion based on 77 fusion into Higgs, decaying into /i + yU~, may be very efficient 
in probing the Higgs sector of 2HDM, even for luminosity 100 pb -1 . It is not clear however if 
these low energy options will come into operation. 
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Future experiments will clarify the status of the general 2HDM with a light neutral Higgs 
particle - the role of HERA in such a study may be very important. 

By contrast, for the MSSM the potential of HERA even with luminosity 1 fb _1 is relatively 
poor, producing only 20-30 events of h and A, if the mass is in the 45-50 GeV range for 
tan /3=30. 
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